State reconstruction of finite dimensional compound systems via local projective 
measurements and one-way classical communication 
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For a finite dimensional discrete bipartite system we find the relation between local projections 
performed by Alice and Bob's post-selected state dependence on the global state submatrices. With 
this result the joint state reconstruction problem for a bipartite system can be solved with strict 
local projections and one-way classical communication. The generalization to multipartite systems 
is straightforward. 
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Quantum state reconstruction relics on the ability to 
measure a complete set of observables and further ma- 
nipulation of numerical data in a way to describe unam- 
biguously a quantum system state. Several approaches 
have been given regarding the reconstruction of sole and 
joint quantum system states OH]- Remarkable experi- 
ments have been realized, employing these reconstruction 
methods (See e.g. Refs. [3j H]). For a discrete Hilbert 
space, quantum state reconstruction is achieved through 
the determination of the complete set of real parameters 
describing the state [5]. For a continuum Hilbert space it 
relies on homodyne tomographic reconstruction [2]. Un- 
fortunately full state reconstruction is highly demanding. 
For a multipartite system composed by N qubits, for ex- 
ample, it is well known that 4 N — 1 different real param- 
eters (the generalized Stokes parameters), and thus an 
irreducible number of 3 N different settings are required. 
Supposing a qubit encoding over light polarization, the 
reconstruction is achieved through 2 x 3 N different mea- 
surements, with 2N photodetectors, most of them re- 
quiring multiple coincidences. Although it is certainly 
not possible to reduce the number of parameters to be 
determined it is important to investigate alternative re- 
construction schemes reducing actual experimental lim- 
itations imposed by multiple coincidence measurements, 
and/or the number of detectors. 

In this paper we propose a scheme to reconstruct the 
state of a bipartite system of arbitrary finite dimension 
using only strict local projections and one-way classical 
communication. This scheme follows naturally from the 
answer to a more general question: if Alice performs pro- 
jective measurements on her state, how this affects the 
description of Bob's state? To answer this question, we 
adopt an operational approach based on the submatrices 
of the global system: after Alice's measurement, Bob's 
state will have some dependence on the submatrices of 
the initial global system. We obtain then a general con- 
dition that an arbitrary set of projections must fulfill 
in order to determine completely the global state. We 
analyze well-known and experimentally accessible pro- 
jections, benefiting from previous derivations of single 
system state reconstruction based on number state pro- 



jectors [3 E], as well as the reconstruction of bipar- 
tite Gaussian continuous variable states employing local 
measurements and classical communication |10j . Focus- 
ing primarily on the simpler and important case in which 
Alice's subsystem is two-dimensional (a qubit), we an- 
alyze the projections arising from her measurements of 
Pauli matrices - which are implementable for photons, 
spin-1/2 particles and two-level atoms [3j SI [5] - obtain- 
ing thus an experimentally feasible protocol for state re- 
construction. We then focus on arbitrary bipartite di- 
mension systems, analyzing a set of projections [5J [7] to 
derive a state reconstruction protocol. Finally, we show 
a straightforward generalization of our findings for mul- 
tipartite systems. The detections performed will not be 
dependent on coincidence verifications, and the number 
of detectors employed can be reduced in some cases. This 
is crucial for the important case of multiqubit state or 
joint correlations reconstruction, and thus for multipar- 
tite entanglement detection [TTj . 

An arbitrary density matrix of a bipartite finite di- 
mensional system shared by two parts, Alice and Bob, 
is recasted in terms of well-suited d B x d B submatrices 
Aij in the computational basis {|0, 0), |0, 1), . . . , |0, d B — 
1), |1,0), . . . , \l,d B -l), • • ■ , \d A ~l, dfl-1)} (with d A and 
d B being the dimensions of Alice and Bob 's subsystems, 
respectively) as 
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The ordering of the computational basis is made in such 
a way as to simplify algebraic manipulations. We now 
derive two important results. 

Lemma 1: Alice's density matrix is given by 
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i.e., Alice's density matrix element (pAjij is the trace of 
global state submatrix A^. 
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Proof. In the computational basis, p = 
Ei,fc=o Y?] B iZo Pijki\ij){kl\. Alice's state is given by the 
partial trace X)^=o 1 b{v\p\v)b- An arbitrary term of this 
sum, b(v\p\v) b , is easily seen to be Ei,fc=o Pivkv\i)(k\ 
and Alice's state is then pa — YlJLa 1 b{v\p\v)b — 
J^t^Q 1 Sf,fc=o PivkvVi (k\. An arbitrary element of p A is 

(PA) m = (p\pa\v) = J2u=o 1 P^w and an arbitrary ele- 
ment of the submatrix A m is (A m ) a p — p^ ari p- Finally, 
the trace is simply TrA m = J2t=o 1 P*u"7" = (P^W ■ 
Proposition 1: Let |t/>) = ES^o 1 Q ™l m ) be a state in 
Alice's subsystem and let Pim — \ip)(ip\ be the projector 
in this state. Then we have that Bob 's state after Alice 's 
projection in state \tf>) is given by 
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that Bob, or Alice, be able to perform local tomography 
of its state. General methods for state estimation such 
as maximal likelihood are standardly applied in imper- 
fect tomography and are inherent to local reconstruction 
methods [TJ [5] . Thus to simplify we assume that Bob is 
able to perform local tomography with Nb copies. 
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Proof. Take P w = Y,tn,nto a m a *n\ m )( n \- 

Then it is easy to check that (P^ <g> Ib)p — 
E^,t=o SjjU 1 a m a* nPnjk i \mj) (kl\. Now, trac- 
ing out Alice's subsystem we get XYa[(P/j <8> Ib)p] — 

Etio 1 a{v\ (EJ^Lo Tfj B i=o a m a* n p n3k i\mj)(kl\j \v) A 

Z^«=o a ma* n Pnjmi\j) We see then that in 

Bob's subsystem TVa[(P|i/,) ® Ib)p] is the same matrix 

as Y^n=a am0l nAnm- We conclude that Bob's state 
after a projection in Alice's subsystem is given by 
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and with the help of Lemma 1 we obtain Eq. ^ , ■ 
(i) General reconstruction projections. We give now 
the general procedures for reconstruction of bipartite 
states. Note that the expression in Proposition 1 can 
be rewritten as 
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In general, the bipartite state reconstruction problem can 
be solved if we use a series of projections on different 
states {tp^) = Ei c 4 l/ ' ) N)' where v — 1,2,...; making 
these projections amounts to constructing the following 
system of equations: 
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Performing suitable projections, so that we can invert 
the above system of equations, allows one to obtain the 
submatrices A mn , and this is the same as to determine 
the state. So the condition to be fulfilled by any set of 
projectors used by Alice is that the above system of equa- 
tions must be invertible. An additional requirement is 



FIG. 1: (Color online) Sketch of the protocol for local re- 
construction of a qubit-qudit joint state, (i) Alice measures 
<t z projecting on |0) or |1) depending on the outcomes ±1. 
(ii) Alice measures a x and a y projecting in |0) ± |1), and 
|0) ± respectively. Alice results are classically communi- 
cated to Bob. 

(ii) 1 Qubit <g> 1 audit. Let us consider first the sim- 
plest case when Alice's state is a qubit and Bob's is 
an arbitrary one (qudit). This is an important example 
since in this case the reconstruction resembles the stan- 
dard procedure for state reconstruction of qubit systems 
[5], apart from the unnecessary coincidence detections. 
Here these are replaced by a communication protocol 
over a classical channel. An arbitrary density matrix 
in the basis {|0, 0), . . . , |0, d B - 1), |1, 0), . . . , |1, d B - 1)} 

is given by p — I +° 01 ), where Aij are rfs x d B 
V Ah fn ) 

submatrices. Also, Alice's density matrix can be written 

PA = i ( 1 + ZA Xa ~ %VA \ in terms of the Stokes 
2 \xa + Wa 1 - z A J 

parameters 1, xa, Vai and za- Alice's Pauli matrices' 
projective measurements amounts to performing the pro- 
jections P| ), P|i>, P|o)±|i) and P| >±i|i>- Equations ^ 
in this case read 



A - 1 + ZA J 0> 
A oo — — ^ — P B > 
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A o + A 11 ±(A 01 + Al 1 ) = (l± XA )p^ ! , (8) 

+ = (l±yA)p l B )±lW . (9) 

The above system of equations for the submatrices is 
easily seen to be overdetermined. This results in a re- 
dundancy, which can be eliminated dealing appropriately 
with the measurement results. Such an issue appears 
naturally in the description of the protocol, which we do 
now. Suppose Alice and Bob share many copies of the 
system whose state they want to reconstruct. Suppose 
also that Bob needs Nb copies of an arbitrary state in 
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his part to perform a good local tomography. Alice starts 
measuring a z on her states, performing projections Piq) 
and Ph\ according to whether outcomes ±1 occur. Alice 
communicates her outcomes to Bob, who separates his 
states in two different subsets, according to these out- 
comes (see Fig.f ). Alice continues the measurements un- 
til both outcomes occurred at least Nb times. Bob then 
by any way performs the local tomographies of p^ B and 

p'g . Note that Alice's outcomes are enough for Bob to 
obtain the value za [H]; by ([7]), Bob obtains the diag- 
onal submatrices Aqq and An. A subtle aspect should 
be noted here. What if Alice measures a z and obtains 
too many outcomes +1 - and hence very few (—1) - so 
that Bob cannot estimate properly p z B ~1 In this case, the 
mean value za = (& z ) is very close to +1: we can assume 
then that An = ~£ A p z B ~ ~ Od B xd B , the null matrix of 
dimension ds X ds- 



The procedure to obtain the off-diagonal submatrix 
.Aoi has a slight modification. Alice measures a x on her 
states, obtaining outcomes ±1, what amounts to respec- 
tive projections P|o)±|i)- Alice communicates her out- 
comes to Bob, who separates his states in two different 
subsets, according to the outcomes. The modification 
here is that Alice will measure until one of the outcomes 
occurred Nb times. When this situation occurs, she can 
stop to measure; Bob then by any way determines the 
density matrix corresponding to this outcome. Using ^ 
and the already determined diagonal submatrices Aqq, 
An Bob obtains Aqi + A' 01 . Once again Alice's out- 
comes are enough for Bob to obtain the value xa\ note 
also that he needs to use only one of the equations Q , 
eliminating half of the redundancy. Repeating this proce- 
dure for u y it should be clear that Bob obtains Aqi — A\ x 
by ([9]) and hence A i; the other half of the redundancy is 
eliminated in this step. Alice and Bob then reconstruct 
the state using only strict local projective measurements 
and one-way classical communication. It is important to 
realize that the parameters xa, Ha and za will be very 
well estimated. For a two-qubit polarization state im- 
plemented in optical systems, Aij are 2x2 submatrices, 
which can be determined by Bob through the same kind 
of measurements employed by Alice. The number of dif- 
ferent experimental settings is 3 2 = 9, as in the usual 
state reconstruction [5] , but detections do not need to be 
in coincidence. As a consequence, that will require only 
2 photodetectors as opposed to the 4 used in the usual 
state tomography [5]. The only requirement for Bob are 
the results ±1 of Alice's measurements, so that he can 
discriminate his states. 

(ii) 1 Qudit £g> 1 qudit. We present now a set of projec- 
tors suited to the case when both subsystems are of arbi- 
trary dimensions. These projectors may be implemented 
in optical systems as demonstrated in [5J [§] or employing 
Stern-Gerlach apparatuses Most importantly, they 
suffice to determine the elements of Alice's density ma- 



trix, which are necessary in the protocol. These projec- 
tors are Pu\ (with j varying from to dA — 1), P|j)±|fc) 
and -P|j)±i|fc) (with j < k and both varying from to 
dA — 1). Explicitly Eq. (|6j) for these projectors writes as 
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(10) 

Aife - {( P A)n + (pA)kk±2Re[(p A ) jk }}p^ ±lk) ,(n) 
fijfc = {{pAjjj + (pA)kk T 2/m[( jOA ), fc ]}^ >±l|fe> (12) 



■A] k ),n 3k = A n +A kk ± 



where A jk = A j:j + A kk ±(A jh 

i{Ajk — Ajifc); with k > j = 0, - ■ ■ ,d A — I. The resulting 
system of equations for the submatrices is once again 
overdetermined. We explain now briefly the protocol, 
which resembles the qubit-qudit protocol in many steps. 

Suppose that Alice and Bob share many copies of the 
state they want to reconstruct and that Bob needs N B 
copies of an arbitrary state in his subsystem to perform 
a good local tomography. To determine a diagonal sub- 
matrix Ajj, Alice starts measuring Pij\ on her states, 
communicating to Bob whether or not (outcomes 1 or 0, 
respectively) the projection actually happened. She re- 
peats this procedure until there happened N B outcomes 
1. Bob then determines by any way the density matrix 
of his states corresponding to outcomes 1. As Alice 's el- 
ement {pA)jj is precisely the mean value (P\j)), Bob can 



get this value with Alice's outcomes and hence, by (10 1 



he determines Ajj. However, if (pA)jj is too small, then 
we can approximate Ajj by the null matrix, as done in 
the qubit-qudit case for the diagonal submatrices. So, we 
can reconstruct all diagonal submatrices. 

To determine a off-diagonal submatrix Aj k (j ^ k), Al- 
ice chooses first to perform one of the projections P\j)±\k) 
on her states. Due to the redundancy in equations (JXTj) , 
only one of them is enough to determine Aj k + A^ k . It 
is important, however, that Alice be able to perform 
the other projection. Alice then perform the projection 
chosen on her states, communicating to Bob whether or 
not (outcomes 1 or 0, respectively) the projection actu- 
ally happened. She repeats this procedure until there 
happened N B outcomes 1. Bob then determines by 
any way the density matrix of his states correspond- 
ing to outcomes 1. Alice's element Re[(pA)j k ] is deter- 
mined if one observes the property Pu\ + ifc\ + Pu\_|j%.\ = 
2 ( P \j)+ P \k))- Th en, we have that (P|j)+|k>)+(Pjj>-|k>) = 
2{pA)jj + 2{pA) kk = gj k , where gj k is a constant value 
which is already determined, since Alice 's diagonal ele- 
ments are known. As Re[(p A )jk] = i p \j)+\k)) ~ ( p \j)-\ty) 
[HI [7], we can take Re[(p A ) jk ] = 2(P| i)+ | fe) ) - g jk or 
Re[{pA)j k ] = gj k — 2(Pj\_H;\). This idea can also be 
used in the situation where the projector chosen has a 
too small mean value, impossibiliting Bob to have N B 
copies on the corresponding subset. We can shift to the 
other projector, which will have a compensating greater 
value. There is, however, a third extreme situation, when 
both mean values are small. But, in this case, (pa)h 



4 



and (pA)kk are necessarily small as well, and the value 
(j>A)jj + {pA)kk ± 2Re[(p A )jk] in (El]) will be too small; 
then we can approximate Ajk + A- k by the null matrix. 



In all situations 
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Finally, Alice chooses to perform one of the projec- 



tions Pi 



\j)±i\k) 



and repeats the above procedure. The 
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Bob to get Im[(pA)jk]- Bob then finds Aj^, — A^ k 
by (12) and hence determines Ajk- As all submatri- 
ces were determined, the state is reconstructed, using 
only local projections and one-way classical communi- 
cation. In this sense, coincidence requirements are un- 
necessary in tomography of compound systems. Also, 
as Alice and Bob acts locally on their respective sys- 
tems, in some cases [12] Alice can eliminate 1 detec- 
tor in her subsystem, due to the completeness relations 

Yfjio 1 P\j) = I, EijUo^V-H^ + P* 



Y,j%=o p \j)+i\k)+P\3)-i\k) = (d A -l)I- where we assume 
j < k. If Bob uses the same set of projectors, he can also 
eliminate 1 detector by the same reason; we get thus an 
economy of 2 detectors. These considerations are even 
more dramatic when considering multipartite settings. 

(in) Multiple Qudits. When dealing with a system 
composed of N subsystems, the state space is composed 
of the tensor product of the individual Hilbert spaces: 
H = Hi ® i?2 <8> • • • ® -Hjv- We can view this system as 
a bipartite one, considering the tensor product of one of 
them with the others 



H= E x 



<g)H N = H 



(13) 



We can then apply the protocol for a bipartite state re- 
construction to this bipartite system, making the pro- 
jections on Hi and sending outcomes to H 2 - After the 
projection, Hi will be left in a pure state and so will not 
affect subsequent operations made on H B (i). Tomogra- 
phy of the various states p^ 



|i>±|fc> |j>±i|fc> t b 



realized in order to perform the task of reconstruction (if 
one uses the projections Pij\, P\j)±\k), P\j)±i\k))- This is 
achieved by repeating the above procedure 



H B (i) — H 2 



Hn — H 



AW 



H 
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(14) 
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After N — 1 repetitions of this procedure we can recon- 
struct the global state density matrix. Coincidence is an 
unnecessary requirement, with an economy of N detec- 
tors for the set of projectors considered. For the case of 
n qubits, one concludes that reconstruction task requires 
n detectors, in contrast to usual In employed, without 
any need of coincidence requirements. 

Conclusions. Quantum state reconstruction is a ex- 
tremely important task in any implementation of quan- 
tum computation or quantum communication protocol. 



It is supposed to determine unambiguously a single or 
joint quantum system state. Unfortunately its implemen- 
tation for composed systems is severely compromised for 
the requirement of joint measurements, which increases 
in number for both increasing the individual Hilbert 
space dimension and/or number of subsystems. We have 
given an alternative procedure for reconstruction based 
on local measurements and classical communication that 
enables bipartite or multipartite quantum systems of ar- 
bitrary Hilbert space dimension to be reconstructed with- 
out the need of joint or coincidence measurements. Be- 
yond the natural importance in quantum computation 
and quantum information, we believe it has important 
implications for the reconstruction of manyparticle sys- 
tem state from simple measurements of local moments, 
such as magnetization for spin systems |llj . The actual 
reconstruction of such a kind of system state through this 
method will be thus a remarkable achievement. 
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This value is given by za = (o"z) = "^ + "~ , where n± is 
the number of ±1 outcomes, respectively. 
This will depend on the experimental implementation of 
the projections. For example, if one uses Stern-Gerlach 
apparatus to generate them [6], the diagonal projectors 
P\j) will be generated by the observable J z . Obviously 
one can then eliminate one detector when measuring the 
diagonal elements, due to the normalization condition on 
the state. However, the off-diagonal projectors are not 
all generated at the same time and thus we cannot make 
such an elimination. 



